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Abstract. We characterize the solution of Navier-Stokes equation as a sto- 
chastic geodesic on the diflteomorphisms group, thus generahzing Arnold's 
description of the Euler flow. 
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1. Introduction 

Laplacian canonical determinism (quoting Laplace himself "We may regard the 
present state of the universe as the effect of its past and the cause of its future... ") 
has been called into question by the studies of H. Poincare and J. Hadamard of 
complex highly unstable trajectories of a system of three bodies moving under 
the effect of gravitation forces. Much later E. Lorenz considered some simplified 
model for atmospheric convection with a small number of degrees of freedom and 
also found chaotic motions, the term chaotic meaning that small changes in the 
initial conditions produce big perturbations in the behaviour of the trajectories of 
motion (sometimes popularly referred to as the "butterfly effect"). 

Dynamics may therefore exibit a chaotic (random?) behaviour even when 
the governing equations are deterministic. Beyond classical dynamics, the most 
convincing illustration of this is, of course, quantum mechanics. 

In Turbulence, modelized by the (deterministic) incompressible Navier-Stokes 
equation 

Ou 

— + u.Vu — vlS.u — Vp, div u = 0, u(0, x) — uo{x) 
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with small viscosity coefficient v > 0, such classical chaotic phenomena really 
occur. On the other hand statistical effects also manifest themselves: they are called 
coherent structures (for example the Jupiter's red spot or Agores anticyclone). We 
refer to [15 for an interesting discussion of these problems. 

There are various possible reasons to look for probabilistic content in Fluid 
Dynamics. The obvious observations that initial conditions are difficult to know 
precisely and that the observable computed in practice have to be statistical quan- 
tities points out to some statistical approach, found, for example, in works of 
Vishik, Fursikov, Foias, etc. ([lOj). 

One can of course add some noise to given deterministic systems and study 
the randomly perturbed ones: then the resulting properties depend on the cho- 
sen perturbation and these studies do not describe the chaotic behaviour of the 
deterministic equations themselves. These approaches have known an impressive 
development in the recent years. 

The intrinsic instability of the hydrodynamical motions, more precisely those 
modeled by the Euler equations (the case where the viscosity is zero) was explained 
geometrically without using Probability theory by V. Arnold ([3j). In his famous 
paper Arnold showed that the solution of Euler equation 



corresponds to the velocity of a flow which is critical for the action functional 



where g denote volume preserving diffeomorphisms of the underlying manifold. Ac- 
cordingly Euler equation can be regarded as the geodesic equation for the metric 
and one derives the instability of the motion strictly from geometrical arguments, 
namely from the negativity of some curvature of the space of diffeomorphisms. 
This is not in contradiction with the formation of the abovementionned coherent 
structures since Arnold describes the so-called Lagrangian flow, namely gt{x) such 
that 



and not the usual Eulerian flow (of velocities). It seems that there is a kind of 
chaotic behaviour of the Lagrangian flow and a more predictable character of the 
Eulerian flow u{t, .). 

We note that this "principle" seems to apply as well to properties like ergod- 
icity. Indeed, there are reasons to believe that the flow of velocities is ergodic (cf. [4] 
or a more recent discussion in fll' ) , invariant measures can be constructed (cf. [2] 
and references therein for the two-dimensional periodic case). On the other hand 
the Lagrangian flow is not ergodic; this has been proved in [7 , for the periodic 
n-dimensional case. 

The statistical effects observed in the evolution of the system are, on the other 
hand, related to the large number of degrees of freedom in presence, a situation 
quite distinct from the dynamical systems studied by Poincare or Lorenz. 



— — h u.Vu — — Vp, div u — 




dtgt{x) 



u{t,gt{x)) 
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Following an idea of Yasue etal ([MlIT^), consisting in interpreting the Lapla- 
cian term of Navier-Stokes equation as origin of the random behaviour of fluid par- 
ticles and the drift as an avereged velocity on the stochastic (Lagrangian) motion, 
one is able to 

a) obtain an intrinsic probabilistic representation of Navier-Stokes solution; 

b) derive it as a solution of a (stochastic) variational principle; 

c) deduce the corresponding (stochastic) equation for geodesies, the associated 
Jacobi equation and therefore consider the stability properties of the motion. 

The methodology thus generalize probabilistically Arnold's approach to the 
case where viscosity is not zero. 



2. Geometry of the diffeomorphisms group 

Let Gy denote the group of volume preserving homeomorphisms on a manifold. 
Throughout this paper we consider this manifold to be T^, the two dimensional 
torus, for simplicity, since we can as well treat the torus in any dimension (cf. the 
remark after Theorem 12. 2|) . On a general manifold the formulation of our results 
may be a more delicate matter. 

The tangent space of the group Gy will be denoted by Gv] it consists of 
vector fields with zero divergence. The norm defines on Gv a canonical Hilbert 
structure. More generally, if C = {g : ^ bijection, g,g~^ G i?*} then for 
s > 2, is a C°° infinite dimensional Hilbert manifold (cf. [9l[l6]). The associated 
tangent space, endowed with the H'^ norm, will be denoted by Q'^ and the tangent 
space for the corresponding volume preserving subgroup Gy by Qy . 

We shall consider Fourier developments, namely, for f G L^, 

f{6) 2 ^ sft/(fc) cos(fc.6l) - S/(fc) sin(fc.6l) 

where / is the Fourier transform of /, 

/(fc) = j^^ exp i{k. 9)d0 

and 

= {k^ (fci,fc2) e ■.ki>0 or fci ^0,k2> 0} 

Notice that / is real iff /(-fc) = f{k). 

We take for orthonormal basis of Qy the following vector fields: 

= n4+T[(fc2COSfc.0)ai - {ki COsk.0)d2)] 

|/c| 

Bl = -^[ik2smk.e)di ~ (fci sinfc.0)92)] 

k <E Z"^ — {(0,0)}, where |fcp ~ k\ + k^ and 9, = together with the constant 
vector fields. 
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When s — we omit the corresponding index in the notation. 
The constants of structure of Qy are 

[Ai,At] = JiL(|fc + + \k- ir'BU) 
[BiBn = + ir'Bi^i -\k- ir'Bi_,) 

[-^k^B'i] = — 2|^|i+i|||s+i (1^ ^1'*^"^^*:+/ ^ 1^ ^ ^l'*^"^^*:-;) 
[d„Bl]^hAl 

One can find a detailed proof of these expressions in [6] for the case where 
s = 0, a computation that extends easily to the general case. For the torus of 
higher dimensions we refer to [7]. 

Define the following functions 

A{\k\\l\\k + l\y+^^^^ ^ ' ' ^" ' 

fls . _ I (lU _ /|2(s+l) _ |l,|2(s+l) I |/|2(s+l)\ 

[k, I] = fci/2 - ^2/1 

The Christoffel symbols for the Levi-Civita connection associated to the H'^ 
metric can be explicitly computed; they are given by 

T%^,At = [k,lKaliBl+i+(3liBl_i) 

ni,Af = [kA{-c^hAUi-Ph^l-i) 

In the case where s — the proof can be found in [5] . 

A similar computation for the space Diff{S^) was done in [T]. 

In order to study the geodesies corresponding to the Euler equation in fact 
only this case is (physically) relevant: the norm (i.e., the energy) has to be 
considered, even if one is lead, for technical reasons, to consider the equations in 
higher order Sobolev spaces, where the norm gives rise to weak Riemannian 
structures (cf. [9].[T6]). 

Concerning the corresponding Ricci curvature (for the H'^ metric) we have 
the following explicit expression 

i 
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and an analogous one for Riccis{Bj) (the curvature is a diagonal operator). 

This result was proved in 7 for the case s = 0. We notice that, in this case, 
the last two terms in each of the Ricci curvature expressions are symmetric, 

and we are there left with only the first terms, which reduce to: 

We observe that the Riccis tensors are given by divergent quantities. 
On the Lie algebra Qy we may define the following process 

with s > 2, where Xk-, yk are independent copies of real Brownian motions. The 
stochastic fiow 

dgit) = {odx{t)){g{t)) 
3(0) = e 

(e being the identity) is a well defined stochastic flow of diffeomorphisms (cf. [I2] 
and [131 15] for exponents s non necessarily greater then 2). Its generator is 

k^O kytO 

Theorem 2.1. The process g{t) coincides with the Brownian motion associated to 
the metric H^. 

The proof follows from the explicit expression of the Christoffel symbols, 
more precisely from the fact that T\s = F^s = 0. 

The Laplace-Beltrami operator on vector fields Z reads, 

k^O k^Q 

where denotes the covariant derivative for the metric. 

Then we have the following fundamental result that enables us to relate the 
classical Laplacian term in Navier-Stokes equations to the stochastic geometry on 
the space of (volume preserving) diffeomorphisms: 
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Theorem 2.2. ([5]) When computed on junctionals of the form F{g){9) = f{g{d)), 
f smooth, we have 



L^F{g) = ^Af 



with c = J2k^o fe|2(3+i) 

The proof can be found in [5] (c.f. also [5]). 

Remark 2.3. We mention that most of the results stated in this paper extend 
without too much trouble (but with a certain amount of work) to the case of a 
torus of dimension higher then two, once we choose a suitable basis of Qv and the 
corresponding ones for the higher order Sobolev spaces. 

Denote by Z"^ a subset of Z"^ such that each equivalence class of the equiva- 
lence relation defined by A; ~ fc' if fc + fc' = has a unique representative in Z"^. We 
should take, for every fc ^ 0, an orthonormal basis of {x £ i?'^ : k.x = 0}, namely 
efc,Q, a — 1, d and define 

Ak^a = ek^aC0s{k.9) 

Bk,a = ek,asin{k.e) 

to obtain (after normalization) the corresponding basis of Qy. 

Furthermore we should assume that J2 k aii^k, aY]"^ is independent of i and 
that X^^, Q,[(efe.Q)*)(efc,Q)-' ) = for i 7^ j in order to obtain the generalization of 
last theorem. 

For the generalization of some of the results in this paper to torus of dimen- 
sions bigger then 2 we refer to [7, and to [5]. 

For regularity reasons, in order to avoid geometrical divergent quantities such 
as the Ricci tensors, we are going from now on to consider a truncated Brownian 
motion {t) defined as 

dx^(i)= {Al dxu{t) + Bl dykit)) 

0<\k\<N 

and the corresponding process (t) with generator L'^'^ . Then, by the same 
arguments of Theorem 2.2. we can show that, when computed on functionals of 
the form F{g){6) — f{g{0)), f smooth, we have 



ith CN = J2 



(kif 



0<|fe|<JV |fcl2(= + i) ■ 



3. A variational principle 

For a time dependent vector field u G L'^{[Q,T];Qv), let g^ be a solution of the 
s.d.e. with values in Gy, 
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dg^(t) = {u{t)dt+J—odx'^{t))g^{t), gu{0) = e 

V Cat 

with generator 

A weak solution of such stochastic differential equation can be shown to exist 
(cf. 0). 

Then, when computed on functionals of the form F{g){9) ~ f{g{9)), f 
smooth, we have 

i^FCg) = uAf + u.Vf. 

Notice that the constant appearing in the diffusion coefficient is chosen to produce, 
after "projection" on pointwise functions, the viscosity parameter v in front of the 
Laplacian. 

For smooth functionals F{g,t) denote 

DF{g{t),t) = lim,^oE^'iFig{t + e),t + e)- ^(^(t), i)) 

For (time dependent) vector fields Z{t,.) e Qv, and using the Ito parallel 
transport (associated to the Levi-Civita connection) Ts^t above the process g{t), 
which is the solution of the stochastic differential equation 

dt[Ts^t]Zt = -(rfe)g(t) o dx^ [Ts^t]Zt 
[Ts^s] = Id 

define 

DZ{git),t) - Hm.^o^E^' {%^(t+,^Z{g{t + e),t + e)- Z{g{t),t))). 
Consider now the following (stochastic) action functional 

^[5] = \e WDgmhdt - ]^E\\Dg{mh 

The Navier-Stokes solution can be characterized as the drift of a critical 
process for this action. More precisely we have. 

Theorem 3.1. ([S]) 

The process g^ (t) is a critical point of the energy functional S if and only if 
u verifies the N. S. equation 

du . „ 

— — h u.Vu = vAu — Vp 
ot 

div u = 0, u{T,0) = ut{0) 
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We precise that critical here means in the space of continuous Gy-valued 
semi-martingales g{t) with ^(0) = e with variations defined as 

{Di)S[e.{v)] = -^^S[e.{v)og{.)]U=o 

where 

et{v),veC\[0,T];g^), v{0,.) = 

deAv) ., , , , 

= v(t, et{v)), eo{v) = e. 



4. Stochastic geodesies equation 

To a variational principle a Euler-Lagrange dynamics can be associated. In this 
framework the corresponding equations are stochastic and correspond to a "regu- 
larized" version of the equation of geodesies for the metric. 

Writting efe for a generic element of the basis of Qy {A/, or B^) and u'^ = , 
define the following truncated approximation of the Riccig tensor 

Rs{^)= E Er:.,e,(rL,eJeO^"-([efc,eJ|eOrf,,n'" 

0<|fe|<Af l,m 

= -2 E E[^'^]'K™<,fc+/?fe,m/?^,fc)«'" 

0<|fe|<Ar m 

^ ^ ' ^ (\k\\m\Y+^ 

0<|fc|<Ar ^' " 

Theorem 4.1. Let s > 0. The vector field u{t, 9) = —u{{T — t),9) is a solution of 
the Navier-Stokes equation if and only if u is a martingale of the process g^ . 
Equivalently, 

DDg^{t) - ^^R^iDg:^) = a.e. 

with 

3^(0) = e, Dg^{T)=u{T,g^{T)) 

or 

Remark 4.2. 

1. The presence of the truncated Ricci term is due to the choice of the Ito sto- 
chastic parallel transport (associated with the Levi-Civita connection). This 

term disappears if wc define the operator D using the so-called "damped" 
parallel transport, which incorporates the effect of the curvature. 



Hydrodynamics, probability and the geometry of the diffeomorphisms group 9 



2. The covariant derivation appearing in the last formula of the Theorem is 
the one associated to the L"^ = H° metric, as it should be, since it is the 
(Euler) non linear term of the equation. On the other hand the stochastic 
perturbations of the Euler equation, responsible for the second order term of 
the Navier-Stokes equation, are related to the metric. 

3. The expression 

C-,N- _ ^rN^u) = vAu 

ZCjv 

may be well defined with each of its separate terms divergent, depending on 
the regularity of the Navier-Stokes velocity function u. 

Proof. The term u'^VekU corresponds, as said before, to the non-linear part of the 
Navier-Stokes equation; so, in this proof, we will be analysing the term coming 
from the stochastic perturbation and, for simplicity, we omit the index s in the 
corresponding Christoffel symbols. 
Writing 



we have 



0<\k\<N 



0<|fe|<Af l,m 
0<|fe|<JV l,m 

We analyse the component in (the one in will give a similar result). 
Wc have to sum the last expression in the indices k and /. The two non zero 
contributions come from 



II = [Bi,B^^f'[Bi,Bt]u'''r. 



Computation of /: 
When I = k + m, 



l^^'^] 17, , |6-n [k,k + m] n,\-'+'^R^ 



2|A:|'^+i|m|^+i' ' 2\k\'+'^\k + m\'+'^ 

+ Im^'B^Ju^'r. 



4|/c|2(«+i)|to|' 
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and, when I = k — m, 



M (|2fc - mr+iB|,_„ + |mr^B^)«^" 



4|fc|2(«+i)|m|«+i 

(notice that B^^ = B^J. 

Concerning //, we have, for Z = A; + m, 



and, for I = k — m. 



[k, m 



2 



Summing up in k the non diagonal terms cancel and 



[k, 



m 



2 



, 2|A;|2(^+i) 

0<|/<;|<Ar fe ' ' 

- Oli.|2(s+l)" " 0|i.|2(s+l)" "•+2^ M - 



2|fc|2(s+i) ^ 2|A;|2(«+i) ^ 2|fc|2(«+i) 

Since cjv = Z^o<|fc|<JV ifcpi+TT = So<|fc|<JV |fc|2(s+i) i 



0<|/s|<iV 

that is, we have the spectral decomposition of the Laplace operator, and the result 
follows. 

The first statement of the theorem, namely the characterization of u{g^) as 
a martingale, comes from Ito formula and the fact that the bounded variation part 
of this process is precisely the expression ^ + u'^Ve^w — C^'^u + ^^R^{u). □ 



5. Final comments 

The probabilistic characterization of the Navier-Stokes solutions described above 
allows to study various of their properties. But one can legitemely ask wether it 
shares any light on the construction of the solutions itself. This would amount 
to ask wether we can solve directly (and in a probabilistic way) the stochastic 
equation 

DDg:^{t) + ^R^{Dg:^) = 
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Stated hke this, we do not have a straightforward answer. Nevertheless one can 
use another probabilistic characterization, in terms of forward-backward stochastic 
equations, based on the observation that the process (t) satisfies not only 

dg^it) = [udt + .[^odxN{t))g^{t), 5^(0) = e 
V Cat 

but also the equation 

dut = -Vp{t,g^{t))dt + XtdxN{t), Dg^{T) = u{T,g^{T)) 

Considering these two equations as a system on {g^ {t),ut, Xt) we have a 
forward-backward system to solve (analytically and numerically) (cf. [8^ and work 
in progress). 
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